Optimal quantum control in nanostructures: Theory and application to generic 

three-level system 
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Coherent carrier control in quantum nanostructures is studied within the framework of optimal 
control. We develop a general solution scheme for the optimization of an external control (e.g., 
lasers pulses), which allows to channel the system's wavefunction between two given states in its 
most efficient way; physically motivated constraints, such as limited laser resources or population 
suppression of certain states, can be accounted for through a general cost functional. Using a generic 
three-level scheme for the quantum system, we demonstrate the applicability of our approach and 
identify the pertinent calculation and convergence parameters. 
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I. INTRODUCTION 

Recent years have witnessed enormous interest in con- 
trolling quantum phenomena in a variety of nanoscale 
systems [1]. Quite generally, such control allows to mod- 
ify the system's wavefunction at will through appropri- 
ate tailoring of external fields, e.g., laser pulses: while in 
quantum optics the primary interest of this wavefunction 
engineering lies on the exploitation of quantum coher- 
ence among a few atomic levels [2], in quantum chem- 
istry optical control of molecular states has even led to 
the demonstration of optically driven chemical reactions 
of complex molecules [3]; furthermore, starting with the 
seminal work of Heberle et al. [4] coherent-carrier control 
in semiconductors and semiconductor nanostructures has 
recently been established as a mature field of research 
on its own. In particular with the advent of semicon- 
ductor quantum dots [5] , sometimes referred to as artifi- 
cial atoms, one now has a system at hand which resem- 
bles many of the atomic properties whilst offering at the 
same time all the flexibility of semiconductor nanostruc- 
tures: experimentally basic quantum-coherence phenom- 
ena such as polarization beating [6] or Rabi-type flopping 
[7] have been demonstrated, whereas theoretically effects 
such as coherent population transfer [8, 9] or entangle- 
ment control [10-12] have been proposed. 

In the last few years this research area has received fur- 
ther impetus from the emerging fields of quantum compu- 
tation and quantum communication [13], aiming at quan- 
tum devices where the wavefunction can be manipulated 
with highest possible precision (quantum gates). This 
high-fidelity quantum-state engineering calls for strate- 
gies which allow an optimal suppression of environment 
losses during gating; self-evidently, such outstanding per- 
formance can only be achieved if the most sophisticated 
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experimental and theoretical techniques for optical con- 
trol of quantum nanostructures are put together. It is 
worth emphasizing that hitherto there exists no clear con- 
sensus of how to optimally tailor the system's control, 
and it appears that each field of research has come up 
with its own strategies: for instance, quantum-optical im- 
plementations in atoms benefit from the long atomic co- 
herence times of meta-stable states, and it usually suffices 
to rely on effective models which can be grasped from the 
solution of simplified level schemes (e.g., adiabatic popu- 
lation transfer in an effective three- level system [14]); in 
contrast, in quantum chemistry the complexity of molec- 
ular states usually does not permit schemes which are 
solely backed from the underlying level schemes, and 
learning algorithms, which receive direct feedback from 
experiment, appear to be the method of choice. Fi- 
nally, coherent control in semiconductor nanostructures 
has hitherto been primarily inspired by quantum-optical 
techniques; however, it is clear that control in future 
quantum devices will require more sophisticated tech- 
niques to account for the enhanced dephasing in the solid 
states; a first step in this direction has been undertaken 
in Refs. [11, 12], where the authors have adopted control 
techniques developed in nuclear-magnetic resonance [15] 
to semiconductor nanostructures. 

In this paper we examine the problem of coherent- 
carrier control in quantum nanostructures within the 
framework of optimal control [16-19]. Here, one starts by 
defining the optimality criteria (the cost functional); in 
general, for a desired quantum-state transition this func- 
tional will depend on the final state, the wish to suppress 
population of certain states during the control process, as 
well as other physically motivated constraints, e.g., lim- 
ited laser resources. The grand strategy then is to min- 
imize this cost functional and to find the optimal time 
dependency of the control fields, which, in turn, govern 
the evolution of quantum states through the underlying 
dynamic equations (i.e., Schrodinger or master equation). 
The calculation of the necessary optimality conditions for 
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this optimization problem results in a system of coupled 
equations, which, for high-dimensional systems, may in- 
volve heavy computations. Yet, the clear-cut advantage 
of this optimization approach is the flexibility to steer 
the control strategies through modification of the cost 
functional, thus rendering this technique ideally for the 
purpose of quantum-state engineering. 

We have organized our paper as follows. Our theoreti- 
cal approach is presented in Sec. II; in Sec. Ill we derive 
the numerical algorithm for the solution of the relevant 
equations. As a first example, in Sec. IV we study the 
optimal control of a generic three-level system. Sec. V 
summarizes our numerical results, and we finally draw 
some conclusions in Sec. VI. 



II. THEORY 

Consider the Schrodinger equation for a n-component 
wave function ip e L 2 (C n , [0, T]): 

? V> = ffV, V>(0) = Vo, (1) 

where the Hamiltonian H = Hq + H t accounts for: Hq, 
the unperturbed system; and H e , the coupling to an ex- 
ternal control field e, where 

\\H €(t) \\<K\\e{t)\\, K>Q (2) 

is supposed to hold for all t G [0, T]; finally, -00 is the 
initial state of the system (h = 1 throughout). Note 
that, strictly speaking, the wavefunction description of 
Eq. (1) is only allowed for an isolated quantum system. 
For the problem of our present concern (control in pres- 
ence of dephasing and relaxation) a more general density- 
matrix description would be required [2] to account for 
the incoherent environment couplings. However, follow- 
ing the procedure outlined in Ref. [20] we observe that 
even in presence of such coupling it is possible to define 
a non-Hermitian Hamiltonian of the form (1), account- 
ing for dephasing and generalized out-scatterings, if, at 
the same time, one introduces a further term which ac- 
counts for generalized in-scatterings. Thus, since we are 
aiming at an optimal control of the coherent time evo- 
lution, i.e., we are seeking for solutions which minimize 
environment losses (see also Sec. IV), we can safely ne- 
glect in-scattering terms, and we are led to Eq. (1), with 
Hq being non-Hermitian. 

In the following we shall consider the problem of de- 
termining the control fields, e E L 2 (C, [0, T]), such that 
Eq. (1) is fulfilled. In so doing we shall be guided by a 
number of further constraints, which, all together, form 
the so-called optimal criteria: firstly, we assert that the 
control sequence brings the system at time T to the de- 
sired state ipd £ C"; secondly, we account for the limited 
laser resources through a minimization of the control field 
strengths; thirdly, we may wish to suppress population 



of intermediate states which suffer strong environment 
losses (see discussion below). More specifically, all these 
constraints are summarized in the cost functional: 

■J(iM) : = \\m-Ml^ + \M\l^ C , [0 ,T]) 
1 ™ 

+ 2 ^Z a MoW 2 mc,[o,T]y ( 3 ) 

j=l 

where the constants 7 > 0, on > are the weighting 
factors, which allow to vary the relative importance of 
the various terms and tpj € L 2 (C, [0, T]) denotes the j-th 
component of tp; the last term of (3) penalizes the occu- 
pation of certain components ipj during the control pro- 
cess. Apparently, further constraints could be added in a 
completely similar fashion. The optimal control problem 
under consideration can now shortly be written as 

rninJ(^>,e) subject to (1). (4) 

We now state that Eq. (4) has a solution. 

Theorem 1. The optimal control problem (4) admits a 
solution $,e) G if^C", [0,T]) x L 2 (C, [0,T]). 

Proof. The above theorem can be verified in a completely 
analogous fashion to Ref. [18] (for details see Appendix 
A and Ref. [19]). ' □ 

To calculate the necessary optimality conditions of first 
order for (4) , we use the method of Lagrange multipliers 
[16] to turn the constrained minimization problem (4) 
into an unconstrained one. For this purpose we define 
the Lagrangian function 

Lty,p, e) = Jty, e) + 5Re (p, - (H + H e ) ^) . 

Here, (4>,ip) = Jq T 4> ' ip*dt, where '*' means complex 
conjugate and the dot '•' denotes the usual vector-scalar 
product in C n . 

Consider the minimization problem: Find tp, p, and e 
such that 

L(4>,P,e) = n inf L(ip,p,e), 
ipex°, P ex t ,eex t 

where X t = L 2 (C", [0,T]) and X° = X t n W> : ip(0) = 
ipo }. Here, necessary conditions for a minimum are ob- 
tained by equating to zero the Frechet derivatives of L 
with respect to the triple (ip, p, e). The following opti- 
mality system is obtained 

iip = (H Q + H e )xP, withV(0)=Vo, (5a) 
ip = {H*+H c )p-q, with ip(T) = tp(T) - ip d , 

(5b) 

e = -MP ■ &~m + * -MP ■ (5c) 
7 oe r 7 dti 

where qj = ctjipj and e = e r + iei. Notice that while 
the state equation (5a) evolves forward in time, the ad- 
joint equation (5b) is marching backwards. The control 
equation (5c) provides the control function. 
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III. NUMERICAL ALGORITHM 



holds: 



In this section we formulate a numerical algorithm that 
solves the optimality system (5a-c) for given initial and 
final configurations xp and ipd, respectively. To solve 
this problem, we apply a gradient-type minimization al- 
gorithm, which first determines a search direction with 
respect to the variable e for both the real and imaginary 
part. Then, a simple step-size procedure is applied that 
guarantees a decrease in the cost functional. For given 
e the search direction is calculated as follows: the initial 
condition for the state equation is given by tpo . Once the 
wave function at t — T is computed, the final condition 
for the adjoint equation is given by ip(T) = ip(T) — ipd- 
Thus, the adjoint variable p can be calculated and the 
gradient of J with respect to e can be computed. 

Assume that the interval [0,T] has been discretized 
into a finite number N step s of subintervals of size St, and 
t rn = (m — 1) St. A discrete state variable at t m is de- 
noted by ip m . To obtain a stepsize that guarantees a uni- 
form decrease in the cost functional we use the Armijo- 
rule with backtracking, [21]. In the sequel we denote by 
J(e) := J(ip(e),e), where ip(e) is the unique solution of 
the state equation for given e. Furthermore we decom- 
pose e into its real and imaginary parts, respectively, i.e., 
e = e r + ie t . The whole optimal control algorithm is then 
specified as follows: 



OPC-algorithm 

1. Initialize e old , < c < 1, v > 1, and (3 > 0. 

2. (a) Solve the state equation iip = (Ho + H t aid) ip 

with ip(0) = -00 (marching forward); obtain 

^new 

(b) Solve the adjoint equation 

ip = (Hq +H c oid)p-q with ip(T) = ip(T)-ijj d 
(marching backwards) : obtain p new . 

(c) Determine a search direction 



where 



-G~ 



i f V r J(e oW ) 
V,J(e oW ) J ' 



,dH 



V r J(e old ) = e° ld - fte[p ■ (^-^)*] new , 

f) M 

V l J(e old ) = ef d - 3?e[p • (— ^)*]™ w 

are the gradients for the real and imaginary 
parts of e, respectively, and G is a positive 
definite matrix. 



3. Determine a step size t such that 



J{e old + t{d r + id % )) < J(e old ) + ct 



V r J(e 
Vj{e oLd ) 



old\ 



d r 



(a) lit = (3 fulfills (6), set (3 := vj3, e new := e otd + 
t(d r + idi), and goto 2, else 

(b) 13 := 13/2, goto 3a. 

Taking in step 2c the matrix G equal to the identity 
matrix leads to the usual gradient method, which can 
happen to converge slowly. Another idea is to use for G 
an approximation to the Hesse matrix of J with respect 
to the real and imaginary parts of e, which leads to quasi- 
Newton methods [21]. 

To determine the evolution of the state variable and 
of the adjoint variable we consider the implicit second- 
order Crank-Nicholson scheme. The advantage of the 
Crank-Nicholson scheme is that it is unconditionally sta- 
ble and it preserves the probability density |^>| in case 
of a coherent time evolution [18]. For completeness, we 
give a brief description of the method. Consider the 
Schrodingcr equation (1). Given the numerical solution 
at time step m, the value of the wave function at the 
next time step, m + 1, is obtained solving the following 
problem for ip m+1 



m+1 



St 

Thus ip m+1 is given by 



where / is the identity in C™. Notice the dependence H m 
from time step due to the presence of the control in H . 
The operator (/ + i^H m+1 ) is a n x n complex matrix 
which is easily invertible (it can be computed analytically 
for small values of n). In case of the adjoint equation 
marching backwards in time the formulae above hold by 
inverting the time direction. 

Accuracy of the solution obtained by integrating in 
time using the Crank-Nicholson scheme (or = 1/2- 
method) is known, and we therefore report only the main 



result. Denote with e r 



il>(t m )-ip m ,m= 1, 



, N, 



steps ; 



the error at each time step between the continuous so- 
lution tp(t) and its numerical approximation tp m . Then, 
assuming e° = 0, there exist constants L and C such that 



|e m | < CSt A 



cxp L 



t n 



1 - LSt/2 



(6) 



where L is a Lipschitz constant and C is proportional to 
sup [0;T] |V>"'|- 

IV. MODEL SYSTEM 

To demonstrate the applicability of our approach, as 
a first representative example in this paper we consider 
the three-level system depicted in Fig. 1, which con- 
sists of: two long-lived states 4>i and tp2, which are 
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i>3 




s 



-0i 



FIG. 1: Prototypical A- type level scheme used in our cal- 
culations: 0i and (j>2 are long-lived states, whereas 03 is a 
short-lived state which is optically coupled to both 0i and 
02 (for details see text); wiggled line indicates relaxation and 
dephasing of state 03. 



energetically separated by some amount 5; a state 03, 
which has a finite lifetime because of environment cou- 
pling (wiggled line). Such A- type configurations have 
a long-standing history in quantum optics and have 
been demonstrated successful in the explanation of many 
coherence-phenomena in atomic systems [2, 14, 22]; more 
recently, similar configurations have received increasing 
interest also in semiconductor quantum dots [8, 20, 23]. 

Within this scheme, the system's time evolution is gov- 
erned by the effective Hamiltonian [14, 20]: 




Time delay (t) 

FIG. 2: Results of our calculations for the simplified model 
system described in the text, and using two laser pulses with 
Gaussian envelopes (with full- width of half maximum r). Be- 
fore the pulse sequence the system is in state 0i; black (white) 
areas correspond to the situation that after the pulse sequence 
the population of 02 is one (zero). Negative (positive) time 
delays correspond to the situation that the ei pulse excites 
the system before (after) the £2 one, and the pulse area is 
defined as dt g(t) (we use fj,i = /X2 = 1 and j a — 0.2r _1 ). 



with 

/ o o mi \ / o o m\ 

Hi = , H 2 = JM2 

\ Mi J \ -i^i -i/j, 2 J 



H n = 




(7) 



where the term — ij a accounts for environment losses 
(e.g., spontaneous photon emissions). Furthermore, the 
coupling to the external field reads: 




(8) 



V. RESULTS 

Assuming that the system is initially prepared in state 
0i, in the following we address the question: what is 
the most efficient way to bring the system from 0i to 
02? Since the direct optical transition between 0i and 
02 is assumed to be forbidden we have to use 03 as an 
auxiliary state; however, intermediate population of 03 
introduces losses through environment coupling. Thus, 
which sequence of laser pulses minimizes the population 
of level 03? 



where \x\ and M2 describe the coupling strengths of states 
0! and 02 to the inter-connecting state 3 (e.g., optical 
dipole matrix elements); note that in Eqs. (7) and (8) we 
have implicitly assumed the usual rotating-wave approx- 
imation [2, 20]. Initial and final states are then given 
by: 







ipd 



iST 



and for the optimality equations (5a) and (5b) we obtain 
e = -^e[P ■ (Him ~ * -^e[P ' 



A. Simplified model system 

To gain insight into the general trends, in the follow- 
ing we shall discuss a somewhat simplified model system 
(results of our complete calculations will be presented 
further below). We assume that the three- level system 
of Fig. 1 is subject to two fields e\(t) = g(t — -y) ex P if* 
and e 2 (i) = g(t + y) cx P^|^ respectively, where g(t) 
denotes a Gaussian envelope with full-width of half max- 
imum t; in addition, we assume that the first pulse e\ 
(centered at time t /2) only affects the 1-3 transition, 
and the second pulse e 2 (centered at time —t a /2) only 
the 2-3 one; such approximation corresponds to the case 
that S is much larger than /lie and M2 e [2, 14]. 
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FIG. 3: Results of our optimal-control calculations for 7 = 
0.001, a — 0: Optimal control (left) and transients of ipi, ip2, 
V>3 (right). 



In Fig. 2 we present results for this simplified model 
system for different time delays t between the two pulses 
and for different pulse areas (as defined in the figure cap- 
tion). As regarding the general trends, we observe from 
Fig. 2 that successful population transfer between states 
4>i and (f>2 can be achieved for both negative and posi- 
tive time delays t a . In the first case, the pulse ei excites 
the system before the e 2 one; here, ei brings the sys- 
tem from (f>i to the auxiliary state and e 2 channels 
the population between <j) 3 and the final state <p2', ap- 
parently, the efficiency of this transfer, which is known 
as the stimulated emission pumping technique [14], be- 
comes maximal when the pulse areas are odd multiples 
of 7r. In contrast, for positive time delays, i.e., when the 
ei pulse excites the system after the £2 one, the popula- 
tion transfer is not achieved through intermediate shelv- 
ing of population; for that reason, the sequence of laser 
pulses is called counter-intuitive, and the whole process 
has been given the name stimulated Raman adiabatic 
passage (STIRAP) [14]. This STIRAP process exploits 
the renormalizations of quantum states in presence of the 
strong laser fields, and population transfer is achieved by 
keeping tp 3 negligible throughout; see Ref. [14] for an ex- 
cellent review. 



B. Optimal control 

We next consider the population transfer for our com- 
plete model system of Eqs. (7) and (8), i.e., e affects both 
the 1-3 and 2-3 transitions, within the framework of opti- 
mal control (see Rcf. [24] for a related optimization anal- 
ysis of the coherent dynamics). As will become apparent 
from our following discussion, we can change the charac- 
teristics of the transfer process between the two limiting 
cases through different penalization of tp 3 , i.e., through 
modification of the weight a = 0:3 in the cost functional 
(3). In the solution of Eqs. (5a-c) we use c = 0.001, 
v = 1.1, and initialize (3 with 0.2 (we checked that our 
results do not depend decisively on these parameters). 
Furthermore, we use \i\ = fi2 = 1, 5 = f0, 7o = 0.01, and 
consider a final time T = 40. To update the matrix G 
in our algorithm we use the Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) formula [21]. Unless otherwise specified, 
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FIG. 4: Same as Fig. 3 but for a = 0.01. 



e = 0.1cxp(i<5t) for the initialization of the algorithm, 
and we stop the iteration if the norm of the gradient is 
less then 10~ 5 . 

Figure 3 shows results for (a) the control field and (b) 
the quantum-state population for a = (i.e., no penal- 
ization of ^3). We observe that the population is chan- 
neled through occupation of the interconnecting 1^3 state, 
in close resemblance to the stimulated emission pumping 
process; indeed, analyzing the Fourier transform of the 
control field, Fig. 3a, we find two strong contributions at 
frequencies — (5/2 and 8/2, where the first one dominates 
at times below 20 and the latter one in the second half 
of the transfer process. 

Further insight into the pertinent calculation parame- 
ters can be obtained from table I which reports the influ- 
ence of 7, Eq. (3), on: the tracking error \ip(T) — ipd\c3, 
i.e., the measure of how accurately the final state is 
reached; the value of the cost functional J; and the num- 
ber iter of required minimization steps. Quite generally, 
we observe that allowing stronger field strengths through 
choice of smaller 7- values results in a more efficient pop- 
ulation transfer, as one would already expect from the 
discussion of the simplified model system: here, the in- 
creased control field allows for faster 1^3 and 3^2 
transitions and, in turn, smaller environment losses due 
to population of $3. We also verified that the on- and 
off-switching of the control fields can be controlled by re- 
placing e in the state equation by £,(t)e, where £(i) is a 
function that smoothly approaches zero at early and late 
times; such additional control might be required to ac- 
count for the limited laser resources in experiment. We 
finally note that the number of required iterations of the 
algorithm significantly increases as 7 decreases, which 
we attribute to the increasing singularity of the optimal- 
control problem as 7 — > + . 

Finally, in Fig. 4 we show the influence of the penal- 



TABLE I: 


Results 


of our c; 


ilculations for different 


values of 7. 
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\i>{T) 




J 


iter 


10- 1 


5.76 


• 10 -3 


3.93 • 10~ 2 


26 


10~ 2 


2.68 


■ 10~ 3 


5.53 • 10" 3 


63 


10 -3 


7.44 


■ 10- 4 


1.11 • 10~ 3 


123 


10- 4 


1.68 


■ 10- 4 


2.44 • 10~ 4 


360 
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ization of the 03-population through finite values of CK3, 
Eq. (3), on the solution of the optimal-control problem, 
which results in a strong reduction of the population- 
transfer time. From the Fourier transform of the control 
field of Fig. 4 we furthermore infer that here the time 
ordering of the two dominant frequency components is 
reversed as compared to the a — case, thus making 
this excitation scenario similar to the STIRAP process. 



VI. CONCLUSIONS 

In conclusion, we have presented a theoretical analysis 
of optimal control of quantum nanostructures. A gen- 
eral solution scheme for the optimization of an external 
control (e.g., lasers pulses) has been developed, which al- 
lows to channel the system's wavefunction between two 
given states in its most efficient way; physically moti- 
vated constraints, such as limited laser resources or pop- 
ulation suppression of certain states, can be accounted for 
through a general cost functional. A computer algorithm 
for the solution of the optimal control problem has been 
derived and analyzed in detail. Finally, we have demon- 
strated the applicability of our approach for a generic 
three-level quantum system, and we have identified the 
pertinent calculation and convergence parameters. 

Apparently, the true strength of optimal control can 
only be appreciated in the study of higher-dimensional 
systems where the control strategies can no longer be 
grasped from simple considerations, which are needed, 
e.g., for the design of quantum gates in future quantum 
registers. There, the high flexibility of our present ap- 
proach, which solely relies on the state equation and a 
general functional accounting for the control constraints, 
renders optimal control as an ideal tool for both theoret- 
ical modeling as well as experimental support. In this 
respect, it will be necessary to develop more efficient nu- 
merical methods for the solution of the bilinear optimal 
control problem of our present concern. Future work will 
also address applications beyond the presently studied 
three-level scheme. 
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APPENDIX A 

In this section we present the proof of Theorem 1 . 

Proof. Step 1: We first calculate an a-priori estimate 
for the solution of the state equation (1). For given 
e G L 2 (C, [0, T}) the equation (1) is a linear ordinary 



differential equation, which therefore has a unique solu- 
tion -0 e iJ^C™, [0, T]). We now write Eq. (1) in integral 
form: 

itp(t) =tp + [ (Hq + H e )ip(s) ds for < t < T. (Al) 
Jo 

Taking the Euclidean norm in C™ on both sides of (Al) 
and using the triangle inequality on the right-hand side 
results in 

\\m\\<\m\+ f\\\H \\ + \\H e \\)\ms)\\ds, 

Jo 

where for the matrices Hq and H e 1 1 • 1 1 denotes the induced 
matrix norm. We can now apply Gronwall's inequality 
[19], and obtain 

U(t)\\ < IhMexpQT (||jr || + ||ir £ ||)ds) forO<t<T. 

Using equation (2) and squaring and integrating the 
above inequality results in 

IM*)||l 2 (CM0,t]) < IkAoll-Kiexp (K 2 + K 3 \\e\\ L 2 (c [0 T]) ) 

with constants K\, K2, K3 > which do not depend 
on e. Furthermore, using the state equation (1) yields 
that for all e e L 2 (C, [0,T]) with | |e| \ L *(c,[o,T]) < K 
the corresponding states tp(e) are bounded in H 1 , i.e., 
IIV'(e)l|H 1 (C'*,[o,T]) < K for some K > 0. 

Step 2: Let (ek)k>i be a minimizing sequence for J, 
i.e. 

Urn J(ip k ,e k ) = inf J(tp,e), 

k^oo {'4 } : u ) satisfies (1) 

where we denote by ipk — VK e fc) the unique solution of 
(1) for given efc. Due to the fact that J(ip(e), e) — > oo as 
||e||z, 2 (C,[o,T]) — ► oo, we get that the sequence (efc)fc>i is 
bounded in L 2 (C, [0, T}). Since the unit ball in a Hilbcrt 
space is weakly compact, there exists a weakly to an 
e e L 2 (C, [0, T]) convergent subsequence, which we again 
denote by (ek)k>i- Step 1 above ensures that the corre- 
sponding sequence (ipk)k>i is bounded in i? 1 (C™, [0, T]); 
thus, again by choosing a proper subsequence 

V> fe i> weakly in H\C n , [0, T]), 

it follows from the Sobolev imbedding theorem [25] that 

ipk -> tp strongly in L 2 (C n , [0, T]) and in C°(C", [0, T]). 

We can now show that (-tjj, e) is a solution of the optimal 
control problem. From the definition of tpk , we have 

iMt) = + / {H + H €k )^ k {s) ds for < t < T. 
Jo 

(A2) 
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Next, we consider the limit in (A2). The weak L 2 - 
convergence of to e implies weak convergence also 
for the complex conjugates, i.e. e£ — ^ e* in L 2 (C, [0, T]). 
Strong convergence of ^ to tp in L 2 (C, [0, T}) allows to 
go to the limit as k — > oo on the right-hand side of (A2). 
Thus we find 



iip(t) = ip 



{Ho + He)j>(8) ds, 



which shows that tp = ip(e), or equivalently that (ip,e) 
fulfills (1). We finally obtain 



1 n 
J($^) =2^( T ) - Vttlc» + Jl|e||| 2(Ci[ o,T]) + | H^Hi 2 (c,[o,T]) 

i=i 

1 1 ™ 

<- lim |Vfe(T) - Vdlc" + X 11 , 111 mf l|efe I li»(c,[o,r|) + o , lim S ^ ' ^ ' l* 2 (C[o,r|) 



inf J{i>,e), 

(tp,e) satisfies (1) 



r 



where we used the lower-semicontinuity of the L 2 -norm. of Eqs. (1) and (3). 
Thus we have proved that (ip, e) is a solution to problem 
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